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In previous papers we have shown that scattering of spacecraft nucleons from dark matter 
gravitationaUy bound to the earth gives a possible explanation of the flyby velocity anoma- 
lies. In addition to flyby velocity changes arising from the average over the scattering cross 
section of the collision-induced nucleon velocity change, there will be spacecraft temperature 
increases arising from the mean squared fluctuation of the collision-induced velocity change. 
We give here a quantitative treatment of this effect, and suggest that careful calorimetry on 
spacecraft traversing the region below 70,000 km where the flyby velocity changes take place 
could verify, or at a minimum place signiflcant constraints, on the dark matter scattering 
model. 



I. INTRODUCTION 

In several recent papers we have explored the possibility that dark matter scattering may be 
responsible for the anomalous geocentric frame orbital energy changes that are observed during 
earth flybys of various spacecraft, as reported by Anderson et al. [l]. Some flybys show energy 
decreases, and others energy increases, with the largest anomalous velocity changes of order 1 part 
in 10^. While the possibility that these anomalies are artifacts of the orbital fitting method used in 
[l| is being actively studied, there is also a chance that they may represent new physics. In we 
explored, through order of magnitude estimates, the possibility that the flyby anomalies result from 
the scattering of spacecraft nucleons from dark matter particles in orbit around the earth, with 
the observed velocity decreases arising from elastic scattering, and the observed velocity increases 
arising from exothermic inelastic scattering, which can impart an energy impulse to a spacecraft 
nucleon. In Q] we constructed a concrete model, based on two populations of dark matter particles, 
one of which scatters on nucleons elastically, and the other of which scatters inelastically, each with 
a shell-like distribution of orbits generated by the precession of a tilted circular orbit around the 
earth's rotation axis. We showed in that this model can give a good fit to the fiyby data, with 
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shell radii in the 30,000-35,000 km range. 

In the present paper we follow up on the brief observation in Q] that if there is a spacecraft 
velocity change as a result of dark matter scattering, there must be a corresponding temperature 
increase arising from fluctuations in the scattering recoil direction. In Sec. II we develop formulas 
for giving a quantitative treatment of this effect. In Sec. Ill we give order of magnitude quick 
estimates, and in Sec. IV we give numerical results based on the model of ^3]- In Sec. V we suggest 
that a thermally shielded, spacecraft based calorimetry experiment could potentially give crucial 
information on the dark matter scattering model for the flyby anomalies. 



II. TEMPERATURE CHANGE ARISING FROM VELOCITY FLUCTUATIONS 

In we considered the velocity change when a spacecraft nucleon of mass mi ~ IGeV and 
initial velocity ui scatters from a dark matter particle of mass m2 and initial velocity U2, into an 
outgoing nucleon of mass m-i and velocity vi, and an outgoing secondary dark matter particle of 
mass 7712 ~ ""^s — ^TTT- and velocity V2 ■ (In the elastic scattering case, one has m2 = m2 and 
Am = 0.) Under the assumption that both initial particles are nonrelativistic, so that << 
c, |'U2| << c, a straightforward calculation shows that the outgoing nucleon velocity is given by 

miUi + m2U2 , . f-,-. 

Vl = ■ -. h WVout • (1) 

mi + 777-2 

Here w > is given^ by taking the square root of 

2 m2m2 1^ ^ .2 , Amm'2 I 2 (mini + 7712772)^ 1 

^ =7 7 w —^(^1-^2) + 1 2c , (2) 

(mi + m2j(mi + m2) mi(mi+m2jL (mi + m2j(mi + m2j-l 

and TJout is a kinematically free unit vector. Denoting by 6 the angle between Vont and the entrance 
channel center of mass nucleon velocity ui — (miui + m2U2) / {fni + rn2) = m2{ui — U2)/ ("7i + m2), 
and assuming that the center of mass scattering amplitude is a function f{6) only of this polar 
angle, the average over scattering angles of the outgoing nucleon velocity is given by 

,^ , miUi+m2U2 , , .,771-^2 ,„s 

(7;i) = ■ h77'(cos6')— — , (3) 

mi + m2 |7ii — 7i2| 



with {cos 9) given by 



,,._ Iodesin9cose\m\^ 



^ The notation t was used in 0] for what we here term w; the change in notation avoids confusion with use of t for 
time. We take ft = 1 in Sec. Ill, while the velocity of light is denoted throughout by c. 



Subtracting ui from Eq. ([3]) gives the formula for the average velocity change used in fl and 0] 

to calculate the flyby velocity change, 

, . m2U2 - m'^ui Q.ui-jh .-s 
iovi) = Vw{cos6)-r-z ^ • (5) 

mi + 7712 1^1 ~ ""21 

However, in addition to contributing to an average change in the outgoing nucleon velocity, 
dark matter scattering will give rise to fluctuations in this velocity, which have a mean square 
magnitude given by 

((^1 - (^^i))') = w'^dvo^, - (cos g) ,!^ = w\l - {cos9f) . (6) 

V \Ul — U2\'' 

This fluctuating velocity leads to an average temperature increase of the nucleon, per single scat- 
tering, of 

m = ^{{n-{n)f) = ^wHi-{cosef) , (7) 

2kB 2kB l_. 

with ks the Boltzmann constant. In analogy with the treatment of the velocity change 5vi in ^j, 

to calculate dT/dt, the time rate of change of temperature of the spacecraft resulting from dark 

matter scatters, one multiplies the temperature change in a single scatter {5T) by the number 

of scatters per unit time. This latter is given by the flux j-ui — U2\, times the scattering cross 

section a, times the dark matter spatial and velocity distribution p{x,U2). Integrating out the 

dark matter velocity, one thus gets for dT/dt at the point x{t) on the spacecraft trajectory with 

velocity ui = dx{t)/dt, 

dT/dt = j d^U2{6T)\ui-U2\(Tp{x,U2) ■ (8) 

Integrating from ti to tf we get for the temperature change resulting from dark matter collisions 
over the corresponding interval of the spacecraft trajectory , 

Tf — Ti = / dt d^U2{6T)\ui - U2\o-p[x,U2) ■ (9) 
Jti J 

In the elastic scattering case, with Am = 0, m2 = m2, the formula of Eq. ([2]) simplifies to 

= ( ^^)\ui - n2? . (10) 

V mi + 7712/ 

In the inelastic case, assuming that Am/m2 and 7772/7772 are both of order unity, Eq. ^ is well 
approximated by 

2 ( 2Am m'2 \ 2 f-,-,. 
w ~ -J- c . (11) 

\ ?77i(r77i + 7772) / 

Since ui and U2 are typically of order 10 km s^^, the temperature change in the inelastic case, per 
unit scattering cross section times angular factors, is larger than that in the elastic case by a factor 
cV|ni|2~109. 
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III. QUICK ESTIMATES 

Before going on to detailed modeling calculations using Eq. ([9]), we first give quick estimates 
using Eqs. (jlOp . and (llip . making the approximations that the dark matter mass rn-2 is much 
smaller than the nucleon mass rn-i, and that {cos 9) in Eq. ([7]) is much smaller than 1. In the 
elastic case, Eq. (4) of tells us that the magnitude of the velocity change in a single collision is 
of order 

|((5i7i)| ~ — l-ui - U2I ■ (12) 
mi 

Taking the ratio of the single collision temperature change to the single collision velocity change, 
and multiplying by the flyby total velocity change ~ lO^^I-uij, we get as an estimate of the total 
temperature change 



in agreement with [2]. 

In the inelastic case, we must take into account the kinematic structure of an exothermic inelastic 
differential cross section, as was done in jsl (but was not correctly done in the estimate given in 
Q]). In the inelastic case, Eq. (5) of [2] tells us that the magnitude of the velocity change in a 
single collision (for (cos 6) > 0) is of order 



K^^;,)!^ V?5^c(cos0) . (14) 
mi 



Writing the inelastic differential cross section near threshold in the form 



- ^i-lfc'fc-l + Ancl(A:')':^COS0 + (15) 



we have 



{cos 9) ~Anelfc7(Ancl^ ^) 



(16) 



with k the entrance channel momentum 



mim2 1^ ^1 1^ ^1 , . 

k = ui — ti2 — m2 ui — U2 , (17) 

mi + 7712 
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and with k' the exit channel momentun, which to leading order in Am is 



k' ~ v^2Amm2C . (18) 

Again taking the ratio of the single collision temperature change to the single collision velocity 
change, and multiplying by the flyby total velocity change ~ 10~^|mi|, we get as an estimate of the 
total temperature change in the inelastic case 



Tf — Tj ~ — — -10 \ui\ — — — — —^ — . (19) 

\{dVi)\ Sinol |ni - M2I "12 

Defining a dimensionless parameter 5inei characterizing the inelastic scattering by 

(m2c)25i„ei , (20) 



^incl / \ 2 I 



-^incl 

the estimate of Eq. ()19p can be rewritten as 

7^/-r.~10-'^^Si,ei^^~0.6xlO^°K5„,eif^l . (21) 
■' 2kB \ui — U2\ \ MeV / 

Thus if 5*11101 is of order unity, the inelastic scattering temperature rise is substantially bigger than 
that from elastic scattering, as already anticipated in the remarks following Eq. (jlip above. For a 
dark matter- nucleon scattering force of range a, one expects i?inoi/^inei ~ o?, giving the estimate 
'S'inei ~ l/{'m'^c^a?') = (A2/a)^, with A2 = l/(m2c) the dark matter Compton wavelength. So 6*11101 
can be substantially less than unity. 

Since {cos 9) ~ kk' /{SinQ\m2C^), the approximation of neglecting {cos 9) compared to 1 is equiv- 
alent to assuming that S'inei >> kk' /{m2cf' = kk'}^. We continue to make this assumption in the 
numerical work of the next section, and note that it is compatible with having a small value of 

•S'inei • 



IV. RESULTS FROM A MODEL FOR THE FLYBY VELOCITY ANOMALIES 

In we formulated a dark matter scattering model for the flyby velocity changes, by assuming 
that inelastic and elastic dark matter scatterers populate shells generated by the precession of 
circular orbits with normals tilted with respect to the earth's rotation axis. By some simple 
substitutions in the computer program used to calculate the velocity change predicted by a given 
set of model parameters, one can calculate the corresponding temperature increase of the spacecraft 
predicted by Eq. ([9]). Combining the elastic and inelastic scattering contributions, the results are 
conveniently written in the form 



Tf -Ti _ I 17120^ 

°K ~ I MeV 



(Zol + Sinel-^inel) , (22) 
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„.h and Z„„ dimensionless numbers givin^«pectivdy the elastic and inelastic contnbutions 

to the temperature increase. From fit 2d of [3], for wfiich the elastic and inelastic shell radii 
are 29,370 km and 34,520 km respectively, we give the fits to the velocity anomalies and the 
corresponding values of Z^i and Zi^^i in Table I. Except for Zei for the NEAR spacecraft, which 
in fit 2d barely intersects the elastic shell, the values of Z^i and Z[^e\ are in accord with the quick 
estimates made in the preceding section. 

V. DISCUSSION AND SUGGESTED EXPERIMENT 

To summarize, if scattering from dark matter gravitationally bound to earth is responsible 
for the flyby velocity anomalies, there must also be spacecraft temperature increases when the 
spacecraft passes through the dark matter region. This suggests two space science investigations. 
The first is to analyze the records of both earth-orbiting satellites and earth-exiting spacecraft, to 
see if there are unexplained temperature anomalies, or at least to place bounds on temperature 
increases. The second is to design a dedicated, compact, thermally shielded calorimetry experiment 
that could be carried as a secondary payload on future space missions, to look for temperature 
increases as the spacecraft traverses the region within 70,000 km of the earth associated with the 
flyby velocity anomalies. 
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TABLE I: Flyby anomaly fit 2d and corresponding Zd and Zinoi values. 





GLL-I 


GLL-II 


NEAR 


Cassini 


Rosetta 


Messenger 


6vA (nim/s) 


3.92 


-4.6 


13.46 


-2 


1.80 


0.02 


a A (mm/s) 


0.3 


1.0 


0.01 


1 


0.03 


0.01 


Svth 


3.90 


-4.6 


13.46 


-2.7 


1.80 


0.020 


Zoi X 10^ 


0.26 


0.85 


0.17 X IQ-^ 


0.89 


0.36 


0.38 


Zinel X 10-4 


0.37 


0.36 


0.74 


0.20 


0.65 


0.58 



The first two lines give the velocity discrepancy (5f a and the corresponding estimated error 
OA reported in The third line gives the theoretical values Svth obtained from the model of 
[3], which has two shells of dark matter, one containing elastic scatterers, the other containing 
inelastic scatterers, gravitationally bound to the earth. The Z^i and Zind values, which give the 
flyby temperature rise when substituted in Eq. (|22|) . are given in the final two lines. (In fit 2d, 
on which this table is based, the radius, Gaussian profile width, and tilt angle of the generating 
circular orbit are respectively 29,370 km, 6678 km, and 0.3902 radian for the elastically scattering 
shell, and 34,520 km, 3030 km, and 1.372 radian for the inelastically scattering shell.) 



